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ABSTRACT 



An approximate solution for the temperature distribution 
in a thin circular disk having variable thermal properties, 
variable thickness, and a variable surface heat transfer co- 
efficient is presented. The method used was to cast the 
governing differential equation into the equivalent varia- 
tional problem and solve using the direct method of Rayleigh 
Ritz . Agreement between the approximate solution and an exact 
solution for the case of a constant thickness disk with con- 
stant properties was excellent. A simple computer program 
which may be used to evaluate the temperature profiles is 
provided . 
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TEMPERATURE DISTRIBUTIONS IN THIN CIRCULAR DISKS WITH VARIABLE 
THICKNESS, THERMAL CONDUCTIVITY, AND HEAT TRANSFER COEFFICIENT 
(A DIRECT VARIATIONAL APPROACH USING RAYLEIGH RITZ APPROXIMATION) 

Applications requiring knowledge of temperature distribu- 
tions in disks are many; probably the most common being turbine 
rotors. Exact solution for the temperature distributions in 
disks with variable properties is at present impossible except 
for a few special cases, however approximate techniques may be 
used which allow solution of these problems. The very powerful 
variational calculus approach used in this report is one of 
these techniques. The excellent agreement between the approxi- 
mate solution and an exact solution for a particular case en- 
courages confidence in the analysis. 

The primary motivation for this analysis was the thermal 
stress problem in turbine rotors which requires knowledge of 
the temperature distribution (for the uncoupled problem). Large 
temperature gradients may be present which result in large prop- 
erty variation. Flexibility in the present solution allows an 
iterative approach to the problem of determining a temperature 
distribution for this complex case. The temperature distribu- 
tion may be determined for an arbitrary distribution of prop- 
ertiesj then the new temperature distribution used to obtain a 
new set of property variation and so forth to the desired de- 
gree of convergence. It is expected that convergence would be 
quite rapid. A second intriguing application for those who 
enjoy treading upon thin ice would be to experimentally deter- 
mine the temperature distribution in the disk (rotating or not) 
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then deduce the required heat transfer coefficient distribution 
to produce that temperature distribution. This approach would 
provide simple instrumentation for rotating disks with complex 
flow conditions existing on the faces and still give reasonable 
guantitative heat transfer information. 

Formulation of the Problem . 

Consider an axisymmetric thin circular disk of outside 
radius R and arbitrary thickness distribution b(r). The only 
requirement which we place is that the disk be thin enough to 
ignore temperature variations in the axial direction. This 
allows us to treat the problem as a normal heat conduction 
problem with a uniformly distributed heat sink. We may write 
the governing equation: 

kv z T - Q s -o 

The heat sink term is expressed as a function of the heat trans- 
fer coefficient on the face of the disk. 

CT-'Qo) 

5 

If the dimensionless variables 9,^ are introduced where T is 

1 1 00 

the film, (or bulk) temperature for the heat transfer coeffi- 
cientj 

6 = t - t 00 > 1 =-^ 

then: 

v 2 e -$8=o 

In cylindrical polar coordinates, for axisymmetrical cases 
and thin disks: 

, 1 41 _ £ h R z a 

n <**1 Tb ~ w 
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Figure 1: Sketch Illustrating Geometry and Boundary Conditions. 



Define : 




The boundary conditions of interest are: 




Another set of boundary conditions would perhaps specify the 
heat flux at*y= 1 however for two reasons this was not done 
in the present study, first the heat transfer coefficient is 
high for flow through a turbine as a result of high velocities 
hence the wall temperature will closely approximate the bulk 
fluid temperature; secondly it is difficult to determine the 
actual heat transfer coefficient for the flow through the turbine 
because of the complex geometry. It should be noted that the 
problem of specified heat flux at rotor tip can be handled in 
a way similar to the present problem. 

The present analysis utilized a general second degree 
polynomial to approximate the distribution of X(.r\) • This was 



This distribution was chosen because of the general ease 
with which a parabola may be fitted to a reasonably smooth set 
of data; certainly property data will show few discontinuities. 
Notice that both constant values and linear distributions are 
subsets of this polynomial. More complex polynomials may be 
used at a correspondent increase in algebraic labor. 

The final differential equation to be solved takes the form: 



written : 
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I attempted to cast this equation into some form which would 
allow direct analytic solution through the various theorems 
available for Bessel type differential equations, none of the 
theorems used were applicable to this form. This is the dif- 
ferential equation which is solved by the approximate methods 
of the variational calculus. 

Formulation of the Variational Problem . 

Allow us to consider only the problem in which the approx- 
imate solution will satisfy the boundary conditions, we may 
then write the functional for a homogeneous differential equation 
of the second order as: 



For treatment of more general cases which allow more freedom 
in the choice of the approximating function, the reader is re- 
ferred to Collatz "Numerical Solutions to Differential Equations". 

The Rayleigh Ritz method may be used to write the solution 
in the general form: 



Several specific polynomials were used in the selection of the 
most applicable one. These polynomials are listed in Table I 
and the resultant temperature distributions are shown by com- 
parison in Figure 2. 

The general approach to solution was to substitute the 
chosen approximation into the functional then carry out the 
indicated integration to obtain a function in terms of the 
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Figure 2: COMPARISON OF TEMPERATURE DISTRIBUTIONS 

USING VARIOUS APPROXIMATING POLYNOMIALS. 
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constants a ■ . These constants were then treated as variables 

1 

of the functional and the extremum obtained in the usual way. 

Exact Solution for a Constant Thickness Disk with Constant 
Properties . 

A check of the general validity of the proposed approxi- 
mate method was made by comparison to the exact solution of a 
constant thickness disk with constant properties. Recall the 
applicable differential equation: 

41® + i-A|_zV\K i Cl 
n an. kt> y 

with boundary conditions: 

for the case ~h. K. _ = const = \ 

te- b 

The simple solution results that: 

© _ 

6 c CJT 1 ) 

I = hyperbolic Bessel function of order zero 

This exact solution was used as a criterion for the accuracy 
of the approximation. As will be seen the approximate technique 
reproduced this exact solution with no distinguishable error. 
Discussion 

The various polynomials used are shown below table I. 

e = ©«,[. i + cj. o-vf )] 
i. e - e^O +■ 

s. 0 = ©oL' + AO-n*) a- ft t'-’? 3 )! 

4. e - BoD +<U>-v\ 1 W>V , y] 

5. e = e„[i +=(. O-'i 4 ) + 



at = 1 
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Notice that (1) is just the first approximation to all the 
other functions. By far the best agreement to the exact sol- 
ution was obtained by number five. It is interesting to note 
that numbers 1 and 2 do not allow the function much freedom 
at the point = 0, while 3 and 4 allows and jS to be arbitrary 
at '7^ - 0, and 5 allows d , and P to be arbitrary at 7^ = 0. 

This seems to indicate that the greater the freedom given the 
coefficients the more rapidly will the solution approach the 
exact. Polynomial number 3 allows for an inversion in the 
temperature profile, a condition which is not physically rea- 
sonable for the particular problem being solved. See Figure 2 
for comparison of the various approximation functions. 

A detailed derivation for the case number 5 is presented 
in the appendix for the case of a parabolic distribution of 
the variable . This was treated in the form: 

- A vf + 4 c 

A, B. and C would be specified by the empirical fit of the ex- 
perimental data for a second degree polynomial for the particu- 
lar problem being studied. It would be expected that the dis- 
tribution would be nearly linear for most practical cases and 
as a result the squared term coefficient would be small. An 
actual turbine might be expected to decrease in width as the 
radius increased, the thermal conductivity may either decrease 
or increase with radius (temperature) depending upon the ma- 
terial, and the heat transfer coefficient would be expected to 
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increase with radius due to the higher relative velocities. 

If the term ^ is examined it may be seen how these properties 
influence the parameter. 

zhR 2 

kb 

% tends to increase with radius and increased heat transfer 
coefficient^ and increases with decreased thermal conductivity 
and thickness. A logical calculation procedure is to assume 
the properties constant for the first approximation then use 
the deduced temperature distribution for further approximations 
The results of several computations for ranges of values of A, 
B, and C are presented in Table II. This table presents the 
coefficients for the expression for 

= A Y\ Z + -VC 

and also presents the coefficients oi , (?> , and tf* for the temper 
ature distribution: 

% 0 = / + dt>-y] 2 ) +/3 +■ ip 0-y 6 ’) 

The computer program also calculates the temperature distribu- 
tion but this vast printout would fill a book and the coeffi- 
cients presented here allow for the same information. 

Example of Table useage . 

The use of the Table is completely straightforward. Sup- 
pose we fit a parabola to \ which as the form: 

\ „ 0 1 -v , \ v\ 4-| 

then we enter the table at the row where A = .01, B = .1, and 
C = 1. and read from the same row d ~ .19829, = .020819, and 
y = .000625. The temperature distribution is then: 

e /& o =14-. 0- A 2 *) f (‘"V) KOOOfe'K>0-^ <# ) 
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- . 6 911 4 L - u , 


- • ** 7 1 •» 6 L* - 


1.000 


- . 19637E+0Q 


-.171022-01 


" • ... J> Zs w -* '3. 


2 . C 0 0 


- . 31770E-00 


- . 43354S-U1 


- . 


5.000 


- . 46203E+00 


• - . 142696+00 


- - 2 6 •# 2 -• E * 0 


10.000 


-.45441E+00 


- . 253616+00 


-•llw.£>9t: + 0G 


.100 


-.24275E-01 


- .62856E-02 




.500 


-.IIQIOE+OO 


-.099316-02 


- . 43631E-Oo 
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• 100 


.001 


1.000 


- . 196386+0 0 


- . 17182E-01 


- . 11153E-02 


.100 


.001 


2.000 


- . 31 769E + 0 0 


- . 43420E-01 


- .37226E-02 


.100 


.001 


5.000 


- . 4620 6E+ 0 0 


- . 14273E+00 


-.26437E-01 


.100 


.001 


10.000 


- . 45439E+00 


- . 25361E + 0 0 


- . 11361E + 0 0 


. 100 


. 010 


. 100 


-.24583E-01 


- . 7H50E-02 


. 85506E-04 


• 100 


.010 


.500 


- . 11029E+00 


-.97772E-02 


- .32025E-03 


.100 


.010 


1 .000 


- . 196476+00 


-.17899E-01 


- . 10213E-02 


.100 


.010 


2.000 


- .31764E+00 


- . 44021E-01 


- . 36993E-02 


.100 


.010 


5.000 


- . 46187E+00 


- . 14304E+0 0 


- . 26557E-01 1 


.100 


.010 


10.000 


- . 454236+00 


- .253636+00 


- . 11386E+00 


.100 


• 100 


.100 


- . 27626E-01 


- . 153196-01 


. 17755E-02 


. 100 


.100 


.500 


- . 112216+00 


- . 17486E-01 


. 99955E-03 


.100 


.100 


1.000 


- . 19732E+00 


-.24998E-01 


-.11600E-03 ' 


.100 


• 100 


2.000 


- .31711E + 00 


- . 49963E-01 


- .34931E-02 


• 100 


• 100 


5.000 


- . 46001E-00 


- . 14618E+00 


-.27770E-01 


.100 


• 100 


10.000 


- .45260E+00 


- .25380E + 00 


- . 11629E+00 


.100 


1.000 


.100 


- . 54991E-01 


- .38956E-01 


. 14416E-01 


.100 


1.000 


.500 


- .129356+00 


' - . 86842E-01 


.104026-01 


• 100 


1-000 


1.000 


- . 20471E-00 


- . 8B997E-01 


.56250E-02 


.100 


1.000 


2.000 


- . 31 1 87 E + 0 0 


- .10363E+00 


- . 40053E-02 


• 100 


1.000 


5.000 


- . 442276+00 


" * . 17407E + 00 


-.41241E-01 


.100 


1.000 


10.000 


- . 43721E+00 


- .25383E+00 


-.14112E+00 


.100 


10.000 


.100 


- . 17606E+00 


- .33418E+00 


- . 97287E-01 


• 100 


10.000 


.500 


- . 20251E+00 


- . 31799E+0Q 


- .11276E+00 


.100 


10.000 


1.000 


- . 23062E+00 


- . 30037E + 0 0 


-.131326+00 


.100 


10-000 


2.000 


-.273556+00 


- .27169E+00 


-.16656E+00 


• 100 


10.000 


5.000 


- , 335966+00 


- .21243E-00 


- . 26490E+00 


.100 


10.000 


10.000 


- . 348276+00 


- . 12640E-00 


-.42659E+00 


1.000 


.000 


.100 


-.23378E-01 


- . 55900E-01 


-.291566-02 


1.000 


.000 


.500 


-.104906+00 


- .53764E-01 


-.589666-02 


1.000 


.000 


1.000 


- .18731E+00 


- .56292E-01 


- .94523E-02 


1.000 


.000 


2.000 


-.30418E+00 


- . 72479E-01 


- . 16872E-01 


1.000 


.000 


5.000 


- . 44677E+00 


- . 14947E+00 


- . 49011E-01 


1.000 


.000 


10.000 


- . 44492E+0 0 


-.23884E+00 


- . 14368E+00 


1.000 


. 001 


.100 


- .23411E-01 


- .55979E-01 


-.290206-02 


1.000 


.001 


.500 


- . 10492E+00 


- , 53838E-0 1 


- .58865E-02 


1.000 


. 001 


1.000 


- .18732E+00 


- .56362E-01 


- . 94462E-C2 


1.000 


.001 


2.000 


- . 30417E+00 


- .72537E-01 


- . 16873E-01 


1.000 


.001 


5.000 


- . 446766+00 


-.14950E+00 


- . 49026E- 01 


1.000 


.001 


10.000 


- . 44490E+00 


- .23884E+00 


- . 143716+00 


1.000 


• 010 


.100 


- . 23708E-01 


- .56693E-01 


-.27800E-02 


1.000 


• 010 


.500 


- .105116+00 


- .54512E-01 


- . 57961E-02 


1.000 


.010 


1.000 


- . 18741E+00 


-.56985E-01 


- . 93916E-02 


1.000 


.010 


2.000 


- . 30413E+00 


- . 730616-01 


- . 168806-01 


1.000 


.010 


5.000 


- . 44658E+00 


- . 14977E+00 


- . 49161E- 01 


1.000 


.010 


10.000 


-.44 4r75E + 00 


- . 23884E+0Q 


-.14395E+00 


1.000 


.100 


.100 


- .266506-01 


- . 63753E-01 


-.159976-02 


1.000 


.100 


.500 


- . 10701E+00 


- . 61179E-01 


- .492806-02 


1.000 


• 100 


1.000 


- . 18831E + 00 


- .63151E-01 


-.88777E-02 


1.000 


.100 


2.000 


- . 303746+00 


- . 78248E-01 


-.16973E-01 


1.000 


.100 


5.000 


- . 44488E + 0 0 


- . 15247E+00 


- .505226-01 


1.000 


.100 


10.000 


- . 443256+00 


- .238816+00 


- . 14644E+00 


1.000 


1.000 


.100 


-.53228E-01 


- . 12699E+00 


.64575E-02 


1.000 


1.000 


.500 


- . 124096+00 


- .12101E+00 


. 39183E-03 


1.000 


1.000 


1.000 


- . 196286+00 


- .118596+00 


- . 66919E- 02 


1.000 


1-000 


2.000 


- . 29975E+00 


- . 12492E + 0 0 


-.20231E-01 


1.000 


1.000 


5.000 


- . 428746+00 


- . 17624E + 0 0 


-.65377E-01 


1.000 


1.000 


10.000 


- . 429016+00 


-.23686E+00 


- . 17175E + 00 


1.000 


10.000 


.100 


- . 17543E+00 


- .32184E+00 


- . 1252 7E + 00 


1.000 


10-000 


.500 


- . 20114E+00 


- .304856+00 


- .141066+00 


1.000 


10.000 


1.000 


- . 22855^+00 


- . 28619E+00 


- . 160026+00 


1.000 


10-000 


2.000 


- . 270666+00 

1 


- .25539E + 00 


- .195996+00 



Cf 

i'l 






I 



1.000 


10.000 


5.000 


- . 33317E + 00 


- . 190206+00 


- . 29612E + 0 0 


1.000 


10.000 


10.000 


- . 34820E + 00 


- . 96837E-01 


- .45947E+0G 


10.000 


.000 


.100 


- . 47405E-01 


- .235116 + 00 


- . 15296E+0 0 


10.000 


.000 


.500 


- . 98951E-01 


-.2103 4-6 + 00 


- .168856+00 


10.000 


.000 


1.000 


- . 153266+00 


- . 18481E + 0 0 


- .137296+00 


10.000 


.000 


2.000 


- . 23557E+00 


- . 14715E + 00 


- .22079E+00 


10.000 


.000 


5.000 


-.35688E+00 


- . 89189E-01 


- . 30881E+00 


10.000 


.000 


10 . 000 


- .39249E+00 


-.31690E-01 


- . 452016+0 0 


10.000 


.001 


.100 


- . 47430E-01 


- . 23511E + 00 


- . 15298E + 0 0 


10.000 


.001 


.500 


- . 98970E-01 


- .21035E + 00 


- . 168876+00 


10.000 


.001 


1 .000 


- . 15328E+00 


-.18482E+00 


- . 187316+00 


10.000 


.001 


2.000 


- .23557E+00 


- .14716E + 00 


- . 220816+00 


10.000 


.001 


5.000 


- . 35688E+00 


- . 89186E-01 


- . 30883E+00 


10.000 


.001 


10.000 


- .392496+00 


-.316736-01 


- . 45204E+00 


10.000 


. 010 


.100 


- . 47654E-01 


- .235186 + 00 


-.153156+00 


10.000 


. 010 


.500 


- . 99140E-01 


- . 210 42E + 0 0 


- . 16904E+0Q 


10.000 


.010 


1.000 


- . 15339E+00 


- .184886+00 


- . 18748E + G0 


10.000 


.010 


2.000 


-.23561E+00 


' - . 14721E+00 


’- .221006 + 00 


10.000 


• 010 


5.000 


-.35682E+00 


- .89158E-01 


- .309076+00 


10.000 


.010 


10.000 


-.39243E+00 


- .31519E-01 


- . 45232E+0 0 


10.000 


• 100 


.100 


-.49878E-01 


- .235766+00 


- .15479E+00 


10.000 


.100 


.50 0 


- . 10083E+00 


- .21103E+00 


- . 17073E + 00 


10.000 


• 100 


1.000 


- . 15455E + 0 0 


- .18548E + 00 


- . 189256+ 0 0 


10.000 


• 100 


2.000 


- .23601E+00 


- .14767E+00 


- . 22291E+00 


10.000 


• 100 


5.000 


- .35629'E + OO 


- .88864E-01 


- .31139E+00 


10.000 


• 100 


10.000 


-.39183E+00 


- .29971E-01 


- .455146+00 


10.000 


1.000 


.100 


-.70669E-01 


- .23929E+00 


- . 17246E+0 0 


10.000 


1.000 


.500 


- . 11669E+00 


• - . 21487E+00 


- . 18881E + 0 0 


10.000 


1.000 


1.000 


- . 16542E+00 


- .189236+00 


- .207926+00 


10.000 


1.000 


2.000 


- . 2398 0 E + 0 0 


- . 15020E + 00 


- .24288E + 00 


10.000 


1.000 


5.000 


-.351296+00 


- .84374E-01 


- . 33521E+00 


10.000 


1.000 


10.000 


- . 38634E + 0 0 


- . 13586E-01 


- .48353E+00 


10.000 


10.000 


.100 


- . 19711E + 0 0 


- .12001E+00 


- . 42299E + 00 


10.000 


10-000 


.500 


- .216846+00 


- . 98633E-01 


- . 440316+00 


10.000 


10.000 


1.000 


- .23847E+00 


- . 7 4023E-01 


- .461236+00 


10.000 


10.000 


2.000 


- . 27338E+00 


- . 30394E-01 


- .501126 + 00 


10.000 


10-000 


5.000 


- .33404E+00 


. 73913E- 01 


- . 61142E + 00 


10.000 


10.000 


10.000 


- .368426+00 


. 219316+00 


- . 78408E+00 



TIME, 0 MINUTES AND 51 SECONDS 
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Appendix A - Details of Computations 



The functional may be written: 



= 





i l-hj (/ +■ & -t- Cy B 




Define V as follows: 

~ j™~ r L + ( I n ) 4- ^ ( 1 - ' | 1 ) + ^ ( I - > j 

The derivative becomes : 

' - — 2 <X \ ^ — '^ 1 ' (?> V'"^ — (c? V V ^ 

Squaring 

CV*)" x 4 oL' 4-- 1 Uc4|3» 4“ (^4 ^ 4 - \Ip |3> ) *4 

+ 3ktf 2 n 10 

Squaring the function. 

^ Z - \ 4 2ol 4* 2, ^3 4 2 4* ol Z t V 4" 2 e2 (3> 4 2 (3 S' 

42ol^ - (lot 4- Zol 2 4- 2ol(2 4 2*4^^ -h 

(c^ 2 - 2(3 ~ 2 - 2 ^ ($ -• 2 ,$*) 4 

(Zds 0 - 2 P - ^ 4* 

(.(V 2, + 2^)^ 4 2 P^iq 10 4- f 2 ^' 2 

Define for ease of algebra. 

3 l 4 Z°L 42.(3 4 2^“ + <2 Z 4(2 Z 4 -^ 4 ZoL§ 4 2.(3^ 4 

s - (a<2 4 2dJ -4 2o<_(3 + 2^.^) 

K 3 = oL z - 2(3 -2(3 Z - 2^(3 -2(3^ 

K 4 = 2^3-2^ -2 ^-2^-2^ 

Ks = (3 2 4 S' 

Kfe = 2 £ 5^ 

K-7 a 



Integrals involving 




may be evaluated : 



18 




The integral involving ( / y y )^ becomes: 

J^7/MVy = *( 2 + §<*? + 3UH 1 ~t-Z(3 z + | ?(3if +3i? z 



Substituting these integrals the functional may be written: 



T(T> = =<; 2 -*-§<*0 4-3oiS> + 2(3 2 + fltf + 3ft 7 - 

■*■( 4 A + 3 b + i c) K, + ( 5 A + 5 b + K 2 . 

+•(& A +-=Vb + ±c) K 3 +Cto^ + ^b +^ c) K 4 

+ Cir A +ttb +i%c) Ks- + (tV A + t? b + ^c) Kb 

4 -C-nsA + rs b +^rc) Kt 

The functional is now differentiated with respect to the 
coefficients 0(, /3 , and TP to obtain the extremal. Define 
for convenience/ 



bi = +36 4 4 : c 

L 3 s grA + 4 b c 



Ls - rfc A •+* rr b +• 
l-7 s reA + -rW b + j^c 



X 

t o 



/ - 

v *> 



L 2 = S’A + 5 ^+ 4 ^ 



fc + 



X 



C 



- To A -*• <3 

L (o a r^A+^B+lkc 



Differentiating with 

35L ~ Z ( 1 +• Li i - Z + L 3 ) ^ 4- z. ( j 4*L\ ~L 2 -L 3 +L 4 ) (3 

4 (3 +£il, -2L 2 , -2L 4 -+£Ls)^ + &L,-2L z v ) 
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defining the coefficients of (%, ft, and $ by M*; m' and M^. 
The derivative may be written: ' 

=• (A[ cL + (3 f Mjl 1 + = o 

differentiation with respect to p yields: 

= (§ +2L, - 2L 2 - 2 L 3 4-ZL4^) o4. + 

+ 2L| - 4\-2 +ZL $){ 3 +- 

(.^ +-2L, -ZL 3 -2L 4 +2U')fl 1 +■ (ZLi - 2 L 3 ) 

Defining M’£, M'^, in a similar way: 

+■ Hz (3 * A7/2T 1 f /♦*/' = <? 

finally differentiation with respect to JP yields: 

5^ = 2. L ( — 2 - 2 L 4 +- 2 L 5 ) oi. 4- 

C 74 

■5 + ZL| -2L3 -2L4 +2.L(e)Q> + 

O ^-2L| -4L 4 4 -ZL 7 ) S> 4- ( 2 b -2/.*) = O 

11 " 11 

Define M i > M 2 5 anc * M 3 a s i ma l ar way 

|| = M,'V + M 2' 1 (3 + /p/frV /w/' = £ 

The system of linear equations to be solved for ^ , /3 , 
and TP becomes: 



MiU 


+ 




+ JP = 


“ AI 4 


/yi/'4 




/^/(3 


+• /*?/ <r ^ 


-A?/' 


A//V 


4 


/<'(3 


* = 


- 



Once the desired values for A, B, and C are known the 
included computer program will evaluate d, ft, and 
according to the above method (see Appendix B). 



APPENDIX B 



Computer Program for Evaluation of Temperature Profiles 

Program TEMP carries out the computations outlined in 
Appendix A. This program is written in Fortran 60 for use on 
a Control Data Corporation 1604 digital computer. 

Coefficients used to describe the parameter in terms of 
a second degree polynomial are given the program as input. 

Vrfl - Arf" + B y\ + c 

Output information consists of the temperature distributions 
as functions of the non-dimensional radius. The coefficients 
of the temperature distribution polynomial are also printed 
out . 

A fixed point variable, N, determines the number of sets 
of coefficients to be computed. Symbols and formats for the 



input data are as follows: 

Variable Fortran Format 

A A 8F10.5 

B B 8F10.5 

C D 8F10.5 

N 110 



A sample printout of the program and typical data are present- 
ed in Table 3. A sample output data sheet is presented in 
Table 2. 
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Table 3 



SAMPLE PROGRAM PRINTOUT 
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t . v. T ORi'lAT ( 3F 
END 



10 . 5 » 3X » J ( 1£11.5,3X) ) 



subroutine parafit (c»3.a) 

this subroutine will fit a parabola to three points 

DIMENSION C(3,3), 3(3) » CM (3,3) , A (3), DETCM(3), V( 3 ) 
DEiC-(C(l»lH<C(2>2)*C(3>3)+C(],2)-*C(2»3)-M’f(3,i) + 
1C(1»3)*C(2>1)#C(3»2)) - '(C(l»l)*C(2,3)*C(3>2) + 

2C ( .1 , 2 ) *C ( 2 » 1 ) #C ( 3 » 3 ) +C ( 1 , 3 ) *C ( 2 > 2 ) *C ( 3 » 1 ) ) 

DO 9 OQ n - i , 3 • • 

4 = < 

DO 905 N - 1 , 3 
DO 905 J- 1,3 
9C 3' CM ( N » J ) = C ( N , J ) 

DO .904 1 = 1*3 

904 CM ( I ,K ) - B ( I ) 

i r m ] V '' * N l = , [ CM 1 1 ’ 1 > * CM < 2 » 2 ) *C,M (3*3) +CM (1,2) *CM (2,3) *CM (3,1) + 
1 CM ( 1 9 j ) *-CM ( 2 * 1 ) * CM ( 3 9 2 ) ) - ( CM ( 1 , 1 ) *CM ( 2 > 3 ) #CM (3,2)+ 

2 CM (1,2) * CM ( 2 , 1 ) * C M (3,3) +CM (1,3) *CM (2,2) *CM (3,1)) 

900 A ( K ) = DE TCM ( K ) /DETC 
RETURN 
END 
END 



• 0 0 u 

• 0 0 0 
. lu 
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10.0 
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